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Abstract The statistics of velocity divergence are studied for an assembly 
of particles that float on a closed turbulent fluid. Under an appropriate def- 
inition of entropy, the two-dimensional lagrangian velocity divergence of a 
particle trajectory represents the local entropy rate S, a random variable in 
time. The statistics of this rate, measured in the lagrangian frame, are col- 
lected over a wide range of values. This permits a severe test of the fluctuation 
relation (FR) over a range that exceeds prior experiments, out to a regime 
beyond which the FR no longer holds. Notably, the probability density func- 
tions (PDF) of the dimensionless divergence a T are strongly non-gaussian. 

Keywords Fluctuation Relation • Turbulence 

PACS 47.27.-i ■ 47.27.ed • 47.52. +j 

This work was supported by the National Science Foundation under grant number 
DMR-0201805. 

M. M. Bandi 

Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 
15260, USA. 

Email: mbandi@lanl.gov. 

Current Address: Center for Nonlinear Studies and Condensed Matter and Ther- 
mal Physics Division, MS K764, Los Alamos National Laboratory, Los Alamos, 
NM 87545, USA. 

J. R. Cressman Jr. 

212 Earth and Engineering Sciences, Pennsylvania State University, University 
Park, PA 16802, USA. 



W. I. Goldburg 

Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 
15260, USA. 



2 



1 Introduction 

It is well known that a system subjected to a dissipative process should lead 
to a strict increase in entropy (S). However when a system consists of a small 
number of particles or modes there is a finite probability that, through ran- 
dom action or correlated behavior, the system's entropy can briefly decrease. 
In microscopic systems, where the number of actual particles is small, or 
for systems with few degrees of freedom, transient fluctuations in work [3lJ 
115] or entropy production |19j have been studied theoretically. Experimental 
investigations into microscopic power fluctuations have been performed in 
resistors [27] for the steady-state fluctuation relation, entropy production in 
colloidal particles for the Integral Transient [39] and Transient Fluctuation 
Theorems [9] , work fluctuations in stretched RNA molecules [3"4[[T3] and me- 
chanical torsional oscillators [18,17] for the Jarzynski Equality and Crooks 
fluctuation theorem. 

On the other hand truly large systems can also display improbable fluc- 
tuations when the dynamics are determined by long-lived temporal and spa- 
tial correlations. Recent theories concentrating on macroscopic systems in 
the nonequilibrium steady state have led to the fluctuation relation (FR) 
of Gallavotti and Cohen [26,16,361. The FR has been tested in a variety 
of experimental systems: Rayleigh-Benard convection [TTll38j . incompressible 
turbulence |10| . and power fluctuations in fluidized granular beds [22] . 

This work is an experimental test of the Gallavotti-Cohen steady state 
fluctuation relation (FR) [26lll2il23il25j . The FR concerns the ratio of the 
probability of entropy increase to that of its decrease. The rarity of entropy- 
decreasing events often severely curtails the range over which the FR can 
be tested [30] . This problem is overcome in the present experiment thereby 
permitting a severe test of the FR over an extended range. 

The experiment is conducted on a macroscopic system, namely a large 
tank of turbulent water. The quantity measured is the local rate of change 
of the entropy S in the lagrangian frame for an assembly of particles that 
float on the surface. When the entropy is defined in an appropriate way, the 
FR makes a firm prediction about the time-dependent probability density 
function of S 1 , averaged over all lagrangian trajectories and over sufficiently 
large r. It is observed that the FR no longer holds at large values of an 
appropriately defined dimensionless "entropy production rate" S. The very 
large positive and negative values of S accessible in these experiments have 
hardly been reachable in prior studies. That may explain why the limitations 
of FR have not been clearly seen before. This saturation of the entropy rate 
has been observed in a simple model system [28] and a hint of it is also seen 
in a study of power fluctuations in a resistor 27J. 

The entropy introduced here is that which appears in the theory of dy- 
namical systems [37] ; it is not defined in terms of the heat input and the 
temperature [T1I35] . Rather S is a measure of the rate of contraction of the 
system in the space of its dynamical variables. In this experiment the entropy 
decrease which takes place in two spatial dimensions, is directly apparent to 
the eye, and is amenable to quantitative study. 
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The measurements are all in the lagrangian frame, so that they depend on 
time only, rather than local space variables. The particles form a compressible 
system, allowing them to coagulate and permitting S to be nonzero. On the 
average this quantity is negative, though it can take on both signs from one 
instant to another. 

Under the driving action of the turbulence beneath them, the floaters 
move erratically, causing their local areal density n(x, y, t) to fluctuate in 
space and time. The resulting motion of the floaters in this experiment is 
perhaps surprising. Even if the initial particle distribution is uniform, it 
does not remain so. Rather, the floaters flee fluid up-wellings and cluster 
into string-like structures around fluid down-wellings. This clustering takes 
place within a fraction of a second [29] . The water on which they move is 
incompressible at all values of z including the surface (z = 0). Hence the 
two-dimensional divergence of the velocity at the surface must be non-zero, 
since V2 • v = d x v x + d y v y = —d z v z ^ 0. As emphasized in Ref. [3], com- 
pressibility does not require that the flow be supersonic. Compressibility also 
appears in subsonic flows in which the system of interest is particles that 
have appreciable inertia [4], such as raindrops in a storm. 



2 Theoretical Background 

The starting point is an entropy S(t) = — J dro n(ro, t) lnn(ro, t) defined in 
real space, where the floater concentration n(ro,t) relates to a probability. 
Here ro = r(t = 0) represents the initial positions of particles. These co- 
ordinates are kept fixed while the system evolves in time to ensure a global 
eulerian average over the entire field of view. The coordinate ro is distin- 
guished from the lagrangian position r(t) which changes at every instant as 
the particle is advected due to underlying turbulence. Using this definition 
of entropy and invoking particle conservation, Falkovich and Fouxon [201121] 
show that the entropy production rate can be expressed as: 

S= dr n(r ,t)cj(r ,t) + / ra(r , t) In n(r , t)v • dS. (1) 

J A JB 

Here the differential area element dro = dxodyo and the dimensionless 
particle density n(ro, t) is the fractional number of floaters in a small area of 
the surface centered about each initial point ro and measured at a particular 
instant t. In the notation used here, w(ro, t) = V2 -v(ro, t) is the local velocity 
divergence spatially averaged about a small area also centered about ro- The 
first term on the right of Eq. |T]) , the area term, represents the average entropy 
rate within the area of observation. It is a global average over the entire field 
of view and hence represents an eulerian measurement. The second term 
(boundary term) takes into account the fact that the floaters can leave the 
field of view A. In the area term in Eq. Jl} w(ro, t) is the local entropy rate 
weighted by a probability measure n(r ,i) that leads to a globally averaged 
(over area A) entropy production rate S. 

It is apparent that if the tracked particles were neutrally buoyant and 
small enough, they follow the fluid into the bulk and form an incompressible 
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system. According to Eq. {!]) S = in the incompressible bulk, in spite 
of energy loss due to viscous dissipation, since w(ro,i) is identically zero 
everywhere. 

In these experiments individual particle trajectories are tracked in the 
lagrangian frame. For each particle trajectory r(t), the lagrangian velocity 
divergence uj(r(t), t), or, equivalently, the local entropy rate, is measured. The 
typical duration that a trajectory spends in the held of view is approximately 
700 ms. 

The FR concerns the probability density of the random variable uj(r(t), t), 
or rather, its value averaged over intervals r. Expressing this local divergence 
in dimensionless units, 

_ 1/; dtu(r(t),t) 

° T ~ TP [2) 

Q' = Urn - [ dt u(r(t),t). (3) 



— ' Jo 

Here J?' is the long-time average of the entropy rate, a quantity that is 
negative, reflecting coagulation of the floaters. The ergodicity of the system 
requires that the value of Q' is independent of any individual lagrangian 
trajectory. The usual sign convention adopted for the entropy rate in FR 
studies is positive. In keeping with this notation, we define the time-averaged 
entropy rate as a T = —o' T and Q = —Q' . The averaging interval r is a crucial 
parameter in the FR; the theory is applicable only if r is much larger than 
the correlation time r c of the velocity divergence uj(r(i),t). If r is too large, 
temporal fluctuations of the lagrangian velocity divergence become too small 
to measure. 

The argument r(t) appears on the right hand side (RHS) of Eq.[2]and Eq. 
[3] above, but not on the left hand side. The lagrangian position changes with 
time as implied by the argument t for the lagrangian position r(t). The time- 
average in RHS of above equations is achieved by averaging over local velocity 
divergence of the fluid at different positions along the particle trajectory at 
subsequent time instants. The time- integral therefore implicitly includes an 
average over the lagrangian spatial coordinates r(i). For this reason, each 
surrogate particle (defined below in section [3]) should be viewed as a local 
probe of the lagrangian entropy rate that changes its position with time. 

The ratio of the probability that a T is positive to the probability that it 
is negative, is, according to the FR, 

In [§±^] = a T Hr. (4) 

This equation is put to an experimental test under steady state conditions. 

According to Eq. (J4|, it is more likely that a T is positive than negative. 
The FR makes no prediction about the form of II (a T ). Nevertheless, the 
present experiment yields this PDF for a rather wide range of tr T and for a 
span of averaging times r (see Fig. [IJ . Note that a T has been defined in such 
a way that its ensemble-averaged value is unity. This merely translates to a 
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statement of ergodicity where the ensemble average of the numerator in RHS 
of Eq. [2] equals the long-time average value Q' in the denominator. 

The systems originally considered for derivation of FR were autonomous, 
and followed dynamics of the form r = f(r), thereby assuring time-reversal 
invariance. In the present experiment, the particles evolve under a non- 
autonomous (random) velocity field r = v(r, t) which is time-dependent. 
Tests performed on the particle tracks showed that time-reversal invariance 
was statistically preserved in this experiment. The PDF of the velocity field 
that provided stochastic forcing for the system of particles was found to be 
symmetric about zero, suggesting an equal likelihood of a particle taking 
positive and negative velocities. It is also worth noting that the FR has also 
been derived for systems subject to stochastic dynamics [32 , 33, 7J , the present 
experiment falls in this category. 

3 Experiment 

The experiment is now described in more detail. The tank, of lateral di- 
mensions 1 m x 1 m, is filled to a height of 30 cm. Steady-state turbulence 
is achieved using a pump (8 hp) that circulates water through a system of 
36 rotating capped jets placed in a horizontal plane at the tank floor. The 
turbulence generates ripples on the surface but their amplitude is small and 
can be neglected due to the fact that the energy injection source is far 
removed from the surface where all experimental observations are made. The 
setup is described in detail in Ref . [H] . 

Hollow glass spheres of mean diameter 50 /xm and specific gravity 0.25 
follow the local surface flow. A laser beam (5.5 W) is passed through a 
cylindrical lens to generate a sheet of light that illuminates the surface. Light 
scattered by the particles is captured by a high-speed camera (Phantom 
v5.0) that records particle positions at 100 frames per second. Steady-state 
measurements are made by constantly seeding floaters from the bottom of 
the tank to compensate for those that are lost from the camera's field of 
view. The particles are small enough and light enough that they can follow 
the local flow. To meet this condition, the stokes number St must be much 
less than unity 0]. In these experiments, St ~ 0.1. 

The camera images are captured every 10 ms over an area spanning 9.3 
cm x 9.3 cm. They are stored in a computer and subsequently fed into a 
particle-tracking program in consecutive pairs to obtain the experimental 
steady-state flow velocity fields. A total of 2040 instantaneous velocity fields 
spanning a duration of 20 s are obtained. There are on the average 25000 
randomly distributed velocity vectors in each velocity field, providing reliable 
spatial resolution down to the sub-pixel level. Parameters that characterize 
turbulence on the surface are listed in Table [TJ 

For the analysis discussed below, the experimentally obtained velocity 
fields are seeded with fictitious (surrogate) particles via computer program- 
ming. At t = Q, the 1024 x 1024 pixel grid of the velocity field is deco- 
rated with a uniform array of surrogate particles 6 pixels apart, providing 
170 2 ~ 3 x 10 4 surrogate particles. The mean distance between real particles 
in the experiment is measured to be 6 pixels. Therefore, the initial uniform 
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separation between the surrogate particles is set at this value. Evolution of 
this initially uniform particle distribution is dictated by the experimental 
velocity fields. The particles are tracked by the program instant by instant 
in the lagrangian frame. Though some particles leave the field of view during 
the evolution time, they are tracked as long as they remain in it. For every 
surrogate particle that leaves A, a new particle is introduced at a random 
spatial point within it, thereby creating a new lagrangian trajectory. Real 
particles are substituted with surrogates in order to achieve particle tracking 
in the lagrangian frame. Lagrangian tracking of real particles is beyond the 
experimental capability of the current scheme. 

The instantaneous velocity divergence field uj(tq, t) is calculated by taking 
component-wise spatial derivatives of the experimental velocity field. As the 
surrogates evolve in time, driven by the underlying velocity fields, the velocity 
divergence along the lagrangian trajectory is obtained at every instant. Thus 
one has a time trace of the lagrangian velocity divergence for all evolving 
surrogates in the flow. Over 330,000 lagrangian trajectories were evolved 
in this experiment, of which velocity divergences were recorded for 80,000 
trajectories and employed in the FR analysis. The experiment runs over 20 
s, corresponding to 37 large eddy turnover times (to = 0.54 s), but most 
lagrangian trajectories leave the field of view within 3 turnover times. The 
long-time average of the entropy rate D is +0.37 Hz. 

4 Results and Discussion 

In this experiment only steady-state measurements of a T are reported. In 
practice, that means discarding the first 200 ms of the lagrangian velocity 
divergence time-trace for each trajectory. This 200 ms interval was previously 
measured [2] to be the transient time period for the global entropy rate 
(area term of Eq. [1]) to reach a steady state. The remaining time-trace of 
lagrangian velocity divergence was broken into windows of duration r. The 
time-averaged value from each window represents the numerator of Eq. (|2j). 
Tracks that disappeared in less than (200 ms + r) were discarded. 

Figure [I] is a plot of ln[77(oy)] and ln[7T(— o T )] vs (oy — l) 2 for four 
values of r in units of r c , the decay time of the velocity divergence correlation 
function, C(r) = (u>(ro, t+r)w(ro, t)). Measurements made in the lagrangian 
frame established that this characteristic time is tc — 0.02 s. If the PDFs 
were gaussian, these plots would be straight lines in this figure. The difference 
m[7T(+oy)] — ln[i7(— ay)] = In[i7(+0y)/i7(— oy)] is positive at all values of oy 
in agreement with Eq. (j4|). 

The parameters characterizing all four sets of PDF measurements are 
listed in Table [2j This table shows that the skewness of all the PDFs is quite 
small. Prior experiments, where the PDF of V2 ■ v was directly measured 
and simulated [14j , show that the velocity divergence takes both positive and 
negative values with almost equal likelihood. In that eulerian study the mean 
value of the divergence of the floaters was found to be almost zero. Although 
the skewness for the PDFs is close to the gaussian value of zero, the non- 
gaussian form of the PDFs is apparent in table [2j where one sees that the 
kurtosis is much greater than the gaussian value of 3. 
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Figure [2] shows the LHS of Eq. (J4]) (open circles) and the RHS of that 
equation (open squares) at the same values of t/t c . The measurements in 
Fig. [2] are at odds with the FR (Eq. Q) in two respects: First, they fall 
above the prediction at the smallest value of r/r c = 5, and lie below it in 
the opposite limit of r/r c = 20, but support the FR at the intermediate 
values, t/t c = 10 and 15. It is possible that at early times (t/t c = 5), the 
averaging interval r is not long enough to satisfy the FR prediction, thereby 
leading to its failure. Slow convergence to the FR prediction for averaging 
intervals r > r c has been observed in other experiments |27j . For long time 
averaging (t/t c = 20), the averaging interval r is already approaching the 
average time (700 ms) the trajectories spend within the field of view. This 
also suggests the possibility of a statistical bias due to preferential weighting 
towards trajectories that spend long times there. 

Second, when a T exceeds 7 or so, the ratio TZ = In [ jj/^ T \ ] saturates 
and may even decrease measurably with increasing o~ T . The convexity of 
lZ(o~ T ) is expected under certain conditions |8J. It is however not clear if 
this explanation can account for the present observations. The saturation of 
lZ(a T ) has also been observed in a recent computer simulation of a relatively 
simple dynamical system [28] . In that study the saturation was attributed to 
exponential tails of the PDF measured for the quantity comparable to a T . 
However in the present experiment, the saturation has not proved amenable 
to satisfactory explanation. 

The present measurements extend to larger values of a r than other ex- 
periments and computer simulations 6 . Within the range (0 < a T < 7), one 
sees that the statistical error is small. As a T increases above 7, the number 
of data points Nb decreases, and the statistical noise is of the order of the 
measured ^J\Nb)- The error bars in Fig. ([2]) are (2 x Nb)- It is worth 
noting that the statistical error in this experiment is negligible out to large 
values of <r T . For a gaussian distribution the probability of observing such 
large fluctuations would be vanishingly small. The fat tails of II(<j T ) for both 
positive and negative fluctuations indicate that large fluctuations are not so 
rare, and indeed aid in testing the FR over an extended range of o~ T . 

It is possible that the PDFs in fig. [T] are perhaps gaussian out to (a T — 
l) 2 ~ 50 (er T ~ 7) thereby trivially upholding the FR within this regime. A 
blown up version of fig. [1] within the range < (a T — l) 2 < 50 is shown in 
fig- [3 There the non-gaussian nature of the fluctuations is clearly apparent, 
and seem fully consistent with the chaotic hypothesis and the FR. To the au- 
thors' knowledge this appears to be the first experiment where non-gaussian 
fluctuations verify the FR over an extended range of dimensionless entropy 
rate (o> ~ 7). 

The FR comes in two flavors, the global fluctuation theorem as was ini- 
tially derived [26] , and the local fluctuation theorem The difference 
between the two FRs lies in whether the quantity of interest is measured lo- 
cally at a fixed spatial point or globally over the entire system. Prior studies 
[5] show there is no a priori reason to expect both local and global FRs to 
simultaneously hold for a given system. 

In the present experiment an analysis of eulerian fluctuations in the area 
term of Eq.Q]is possible in principle, and would constitute a test of the global 
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FR. However eulerian measurements capture only large scale flow structures 
under the random-sweeping hypothesis [40] . Although the velocity divergence 
u>(ro, t) is still a local spatial measurement, the concentration n(ro, t) is a pas- 
sive scalar capable of capturing only the large scale flow information within 
the eulerian frame. The eulerian correlation time of the area term (Eq. [T]) 
therefore would become comparable to the large-eddy turnover time (tq ~ 0.5 
s) of the underlying turbulent flow. Since the camera recording is capable of 
spanning only 37 turnover times, a test of the global FR is not feasible for 
technical reasons. Such an eulerian test of the global FR would be possible 
with a direct numerical simulation of the free-surface flow, and it would be 
interesting to learn if the system of floaters would uphold the global FR. 

On the other hand measurements made in the lagrangian frame are ca- 
pable of capturing local flow structures. Accordingly the correlation time of 
lagrangian velocity divergence in this experiment is short enough to make 
an analysis of the FR tractable in the lagrangian frame. In prior tests of the 
local FR [TTlfTU] time-traces of the quantity of interest were collected at a 
fixed spatial point (eulerian measurements). In this experiment, that local 
probe is a surrogate particle which does not remain stationary in space. In- 
stead it changes position from one instant to the other along its lagrangian 
trajectory as it is advected by the local turbulent flow. Within a given av- 
eraging time-window r, the surrogate particle travels a certain distance in 
space, collecting information of the local velocity divergence as it moves. 



5 Summary 

An experiment is reported where the local entropy rate oy, a random vari- 
able, was measured in the lagrangian frame for a system consisting of an 
assembly of particles floating on a turbulent tank of water. The floaters 
form a compressible system, hence exhibit a nonzero value of the velocity 
divergence, which is equal to the entropy rate as defined in the work of 
Falkovich and Fouxon [2"0"ll21| . The definition of a T used here is not that 
of equilibrium thermodynamics but rather is closely related to that used in 
studies of chaotic systems [P6]l37] . The calculations and measurements refer 
to the two-dimensional coordinate space in which the floaters move, not a 
four-dimensional phase space or the full phase space of the large number of 
floaters. 

The experimental steady-state results agree with the fluctuation relation 
over a limited range that the dimensionless entropy rate spans, but a range 
that exceeds that achieved in prior studies. The PDFs themselves are notable 
in several respects: (1) they are only slightly skewed, (2) their mean values are 
small compared to their standard deviations from the mean, and (3) their 
shapes are different for positive and negative values of the dimensionless 
entropy rate a T defined in Eq. $2$. As the coagulation of the floaters implies, 
positive values of the PDF are larger than the negative values at each entropy 
rate. 
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Fig. 1 In [II(a T )] vs. (a T — l) 2 plotted for four values of r in dimensionless units 
t/t c (t c = 20 ms) (a)r/r c = 5, (b) t/t c = 10, (c) t/t c = 15, and (d) t/t c = 20, for 
positive (open circles) and negative (open squares) values of the entropy current 
<j T . For a gaussian the decay is linear. The open circles correspond to coagulation. 



Table 1 Parameters of the turbulence measured at the surface. 



Parameter 



Expression 



Measured value 



Taylor microscale A (cm) 
Rex 

Integral Scale lo (cm) 
Dissipation Rate £diss (cm 2 /s 3 ) 
Kolmogorov Scale r\ (cm) 
Large Eddy Turnover Time to (s) 

RMS Velocity v rm s (cm/s) 



({dvx/dx)*) 



J <(«||M)"> 

„_ (1)1/4 



'/ 

TO = 



= a/<^>-(^ii) 5 



0.3 

93 

1.2 

10.7 

0.02 

0.54 

2.6 



Table 2 Parameters of the PDFs in Fig. \T\ 



Statistics 


T = 5t c 


T = 10t c 


t = 15t c 


t = 20t c 


Mean 


1.03 


0.96 


1.01 


0.87 


Median 


0.94 


0.85 


0.88 


0.75 


Std. Deviation 


9.99 


8.45 


5.49 


5.74 


Skewness 


-0.01 


0.16 


0.18 


0.56 


Kurtosis 


12.84 


16.18 


10.00 


12.98 
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Fig. 2 In [n(+a T )/II(— o>)] (open circles) plotted against o> for the four dimen- 
sionless integration times (a) r/r c = 5, (b) r/r c = 10, (c) r/r c = 15, and (d) r/r c 
= 20. The solid line through the open squares is the RHS of Eq. Q 
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